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ABSTRACT 

Winds from star-forming galaxies provide the most promising explanation for the enrichment of the intergalactic 
medium with heavy elements. Theoretical and observational arguments indicate that the pollution may have 
occurred at z > 6; however, direct observational tests of such a scenario are needed. We model starburst winds in 
the high-redshift universe and hnd that the fraction of space hlled by enriched material varies strongly with the 
assumed star formation efficiency /* and the fraction of supernova energy powering each wind, /esc- We show 
that metals carried by these winds can be seen in absorption against bright background sources, such as quasars 
or gamma-ray bursts, in narrow lines with characteristic equivalent widths ~ 0.5 ^ ^ 5 A. We argue that a 

substantial fraction of the metals in high-redshift winds are likely to reside in low ionization states (C II, O I, Si II, 
and Ee II), but higher ionization states (C IV and Si IV) could also provide useful probes of the winds. The number 
of such lines can constrain both /* and /esc- Statistics of metal absorption lines can also be used to identify whether 
H2 is an efficient coolant in the early universe and to study the initial mass function of stars at high redshifts. 

Subject headings: cosmology: theory - intergalactic medium - galaxies: high-redshift 


1. INTRODUCTION 

It is now obvious that mechanical feedback from galaxies and 
quasars has disturbed a significant fraction of the intergalactic 
medium (IGM). Types of feedback include relativistic quasar 
jets, winds from quasars, and winds from star-forming galax¬ 
ies. These mechanisms have had a variety of consequences, 
including IGM heating (e.g., Voit et al. 2002), the suppression 
of star formation in the host galaxy (e.g., Springel & Hernquist 
2002) and in nearby galaxies (e.g., Scannapieco, Thacker, & 
Davis 2001), and the generation of intergalactic magnetic helds 
(Kronberg, Lesch, & Hopp 1999; Eurlanetto & Loeb 2001). 
Perhaps most importantly, feedback is responsible for enrich¬ 
ing the IGM with heavy elements. Observations of the Lya 
forest at z ^ 3 show that absorption systems with column den¬ 
sities Ahi > 10'“* ^ cm“^ (corresponding to overdensities 5>5\ 
Schaye 2001) have metallicities Z ^ 10“^-10“^ Zq (Tytler et 
al. 1995; Cowie et al. 1995; Songaila & Cowie 1996; Ellison 
et al. 2000). Absorbers with Ahi < 10'"*'^ cm“^ require stack¬ 
ing or pixel-by-pixel analysis because of the weakness of the 
corresponding metal features, but recent studies suggest sim¬ 
ilar metallicities down to systems near the mean density of 
the IGM (Cowie & Songaila 1998; Ellison et al. 1999, 2000; 
Aguirre, Schaye, & Theuns 2002). These observations require 
that > 5% of space be enriched with heavy elements by z 3 
(Madau, Eerrara, & Rees 2001). The era at which this pollu¬ 
tion occurred is uncertain. Songaila (2001) found no evidence 
for evolution in the total mass-weighted abundance of CIV and 
Si IV between z 1.5-5 (see also Qian et al. 2002). However, 
such mass-weighted measurements constrain only the metals in 
the most massive systems and tell us little about the total vol¬ 
ume of enriched space. Furthermore, conclusions about the to¬ 
tal metallicity based on single ionization states are fraught with 
uncertainty. 

The history of metal enrichment has several important conse¬ 
quences for structure formation. Most directly, its extent mea¬ 


sures the volume over which galaxies have influenced their sur¬ 
roundings hydrodynamically, which in turn dehnes their feed¬ 
back on the Lya forest. Additionally, an early phase of metal in¬ 
jection may qualitatively change the characteristic mass of star 
formation (Bromm et al. 2001; Schneider et al. 2002; Mackey, 
Bromm, & Hernquist 2002). The transition between zero- 
metallicity star formation and “normal” star formation has im¬ 
portant implications, e.g., for the expected redshift distribution 
of gamma-ray bursts (GRBs; Bromm & Loeb 2002), for reion¬ 
ization (Wyithe & Loeb 2002d; Cen 2002), and for the chemical 
abundance patterns of low-metallicity stars (Qian & Wasser- 
burg 2001, 2002). 

Stars produced metals inside galaxies, and one must invoke 
some type of outflow to transport them throughout the IGM. 
Simple dynamical removal cannot pollute the large regions we 
observe (Aguirre et al. 2001c; but see Gnedin 1998). Radia¬ 
tion pressure-driven dust outflows may suffice (Aguirre et al. 
2001a), although their effectiveness depends on many untested 
assumptions about the outflows, such as the geometry of the 
magnetic field. Quasar winds are unlikely by themselves to 
entrain enough metals to account for the observed metallici¬ 
ties, although they can help to disperse the metals over large 
volumes. Winds from star-forming galaxies are therefore the 
most likely mechanism. Correlated supernovae in starbursts 
can power “superwinds” with velocities ~ 100-1000 km s“' 
that reach far outside the host galaxy. Such winds have been 
studied in many local starburst galaxies (Lehnert & Heckman 
1996; Martin 1999) as well as in Lyman-break galaxies (LBGs) 
at z ~ 3 (Pettini et al. 2001,2002; Shapley et al. 2002). Because 
their velocities exceed the escape speed from many of the hosts, 
we expect a subset of these outflows to penetrate into the gen¬ 
eral IGM. Escape is most likely from dwarf galaxies, because 
the amount of energy available to power a wind is proportional 
to the halo mass M/, while the gravitational potential confining 
the wind increases with M|. Thus, starburst winds can naturally 


* Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138 
^ Institute for Advanced Study, Princeton, NJ 08540 
^ Guggenheim Feiiow 


1 


2 


account for a high-redshift episode of metal enrichment, when 
the characteristic mass of galaxies is small. 

A great deal of theoretical work has gone into studying 
the galactic wind pollution mechanism. Ideally, the enrich¬ 
ment process should be modeled with cosmological simula¬ 
tions. To this end, several numerical simulations have incor¬ 
porated galactic winds (e.g., Scannapieco, Thacker, & Davis 
2001; Springel & Hernquist 2002; Theuns et al. 2002). Un¬ 
fortunately, we do not yet fully understand the processes that 
generate winds even in nearby, well-observed starburst galax¬ 
ies, although simulations are improving (Mac Low & Ferrara 
1999; Mori, Ferrara, & Madau 2002). We must therefore insert 
simple parameterizations of the wind into the simulation. In 
order to understand the significance of this choice, most work 
to date has proceeded using analytic or semi-analytic methods 
(Nath & Trentham 1997; Ferrara, Pettini, & Shchekinov 2000; 
Madau, Ferrara, & Rees 2001; Scannapieco, Ferrara, & Madau 
2002). Such models have the advantage of allowing much more 
flexible parameter studies. Because the relevant parameters are 
poorly constrained at present, the IGM filling factor (i.e., the 
fraction of space filled by enriched material) estimated by these 
studies varies widely (between ~ lO"”* and ^ 0.5). Aguirre 
et al. (2001c) made a useful compromise between these two 
approaches by inserting winds into already completed cosmo¬ 
logical simulations (see also Aguirre et al. 2001a,b). This al¬ 
lowed a broader parameter study than complete simulations and 
retained many of their advantages (such as the clustering of 
sources), while sacrificing a detailed treatment of the dynam¬ 
ical effects of the winds. 

The crucial unknown parameters in galactic wind models fall 
into three broad categories. First, the characteristic redshifts 
and masses of the wind hosts are critical for setting up simula¬ 
tions with finite resolution and box size. Second, the unknown 
star formation efficiency at high redshifts has dramatic effects 
on both the sizes of individual wind bubbles (because it deter¬ 
mines the energy input into the winds) and on the total amount 
of metals produced. Finally, a variety of questions about the 
structure of the winds remain unsettled. We do not know what 
fraction of the supernova energy is lost to radiation before en¬ 
tering the wind, what fraction of the galaxy mass is entrained 
by the wind, or the geometry of the winds. Studies of nearby 
starbursts can help to answer these questions; however, winds 
at high redshifts may have dramatically different characteris¬ 
tics. For example, local superwinds are bipolar, with the wind 
flowing perpendicular to the galaxy disk. If star formation pre¬ 
cedes disk formation in the halo, we would expect high-redshift 
outflows to be more isotropic (Madau, Ferrara, & Rees 2001). 

Ultimately, we can only answer these questions through ob¬ 
servations of the pollution mechanism. One promising tech¬ 
nique is to study absorption lines caused by metals in the 
winds. In both local starbursts (Phillips 1993; Lequeux et al. 
1995; Heckman & Leitherer 1997; Heckman et al. 2000) and 
LBGs (Pettini et al. 2001, 2002; Shapley et al. 2002) absorp¬ 
tion features can be observed through moderate resolution spec¬ 
troscopy of the host galaxy. Unfortunately, the hosts of high- 
redshift (z > 6) starbursts will be too faint to serve as useful 
background lights. We therefore must use other bright back¬ 
ground sources, such as quasars and GRBs, against which ab¬ 
sorption can be seen. Such background sources have the dis¬ 
tinct advantage of allowing us to probe “relic” winds whose 
host galaxies no longer actively form stars. Relic wind bubbles 
are older and larger than the young winds in starburst galax¬ 


ies; they are also more common than active winds because only 
a fraction of galaxies host starbursts at any time. Thus relic 
bubbles fill more space than active winds, and a careful study 
of their properties is crucial. Similar observations have been 
attempted for local starbursts (Norman et al. 1996) and an ex¬ 
tension of this technique has been used to study the effects of 
LBGs on the IGM (Adelberger et al. 2002). Oh (2002) has also 
proposed using metal absorption lines to study the reionization 
process, provided that metals have been injected relatively uni¬ 
formly in dense regions by z ~ 6. 

The neutral fraction of IGM hydrogen is relatively large 
(> 1%) before the completion of reionization at ^ 6 (Becker 
et al. 2001; Fan et al. 2002). The resulting absorption blueward 
of the Lya transition (Gunn & Peterson 1965) will render unob¬ 
servable any metal lines with rest-frame wavelength Am in the 
range (1 -l-Zr)/(l+z) < ^m/^hya < 1, where ALya = 1216 A and 
z is the redshift of the absorbing system. A transition below 
the left-hand limit will be observable provided that it can be 
separated from the Lya forest. A transition with Am > ALya is 
unaffected by neutral hydrogen inside a finite redshift interval 
Zmin < z < Zs, where Zs is the redshift of the background source 
and 

1+Zmi„ = (l+Z.)^. (1) 

Within this range, metal absorption lines will remain visible. A 
collection of sources with a range of Zs should therefore allow 
us to probe winds throughout the high-z IGM. 

In this paper, we use a simple semi-analytic model for winds 
from star-forming galaxies to study how observations of metal 
absorption by wind bubbles in the era before reionization can 
constrain the metal enrichment process. We compute the halo 
star formation histories and follow the expansion of their winds 
into the IGM. We show that metals carried by the winds produce 
observable absorption features when seen against bright back¬ 
ground sources, and that with current and future instruments, 
observations of these lines can constrain the star formation effi¬ 
ciency in the early universe and can distinguish between differ¬ 
ent halo cooling scenarios and stellar initial mass functions. 

We discuss our supernova wind model in §2 and metals in 
the high-redshift IGM in §3. We present our predictions for 
the statistics of metal absorbers in §4 and discuss our results in 
§5. We assume the standard ACDM cosmological parameters 
of Dm = 0.3, Ha = 0.7, O* = 0.05, erg = 0.9, n = 1, and a Hub¬ 
ble constant Hq = lOOhkm s“^ Mpc“^ We set h = 0.7 in our 
numerical calculations. 

2. SUPERNOVA WIND MODEL 

In this section we describe our model for superwinds in de¬ 
tail. We discuss the abundance of wind sources and their star 
formation histories in §2.1. We fix the mechanical power driv¬ 
ing the winds in §2.2, and we describe our model for their ex¬ 
pansion into the IGM in §2.3. 

2.1. Halo Abundances and Star Formation Histories 

At any redshift, the halo mass function dn/dM determines 
the number of wind sources; each star-forming halo has the 
potential to drive such a wind. We use the Press & Schechter 
(1974) formalism to calculate the mass function, with the mod¬ 
ifications suggested by Sheth & Tormen (1999) to better match 
numerical simulations (Jenkins et al. 2001). 

A halo can only form stars if its mass exceeds a threshold 

Afniin(z) = max[Mfii(z),Mcooi(z)] ■ (2) 
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The first condition requires that the halo mass exceeds the “fil¬ 
ter mass” Mfii, or the effective Jeans mass in the evolving IGM 
(Gnedin & Hui 1998). The second condition requires halo gas 
to cool efficiently after virialization. In the low metallicity ha¬ 
los we study, line transitions of atomic and molecular hydrogen 
are the only effective coolants. Atomic cooling is efficient in 
halos with Tyir > lO'* K, where the virial temperature Tyir is de¬ 
fined in equation (26) of Barkana & Loeb (2001). In most of 
our work, we will fix Mcooi by requiring that atomic cooling can 
proceed. The vibrational lines of H 2 make it an efficient coolant 
in halos with Tyir > 400 K. It is currently unknown when (and 
even whether) H 2 is destroyed by photons below the ionization 
threshold of atomic hydrogen (Haiman, Rees, & Loeb 1997; 
Haiman, Abel, & Rees 2000), or whether radiative feedback 
from early generations of stars enhances the formation of H 2 
(Ricotti, Gnedin, & Shull 2002). We therefore also show re¬ 
sults in which H 2 cooling is permitted. 

The power driving a wind is ultimately determined by the star 
formation history of the host halo, which is usually associated 
with the halo merger history. We use a Monte Carlo imple¬ 
mentation of the extended Press-Schechter formalism (Lacey 
& Cole 1993) to choose histories for each halo. Given a halo of 
mass Ml, at redshift Zf, this algorithm determines the probabil¬ 
ity distribution of progenitor masses (Mi, M/,-Mi) that merged 
to form the object in a small redshift interval dz around Zf (we 
choose dz such that the probability of a major merger is < 1%). 
For each halo, we construct the merger history by first choosing 
Ml from this distribution. We then choose the larger of the two 
progenitor masses (i.e.. Mi if Mi > M/,/2 and M/,-Mi other¬ 
wise) and repeat the process about a redshift Zf + dz. We con¬ 
tinue this procedure until the progenitor mass is smaller than 
Mmin at a redshift Zform- 

In order to simplify our calculations, we divide the merger 
history into discrete events in which the halo mass has grown 
by > 50%. (We find that our results are insensitive to the value 
of this threshold.) We assume that in each such event a fraction 
/* of the accreted baryonic mass Mg is transformed into stars 
on a timescale fsf. We assume that each merger occurs over the 
free-fall time of the host halo fff, and we normally set fsf = %. In 
the rare occasions in which another merger event begins within 
fff of the first, we set fsf equal to the time between these events. 
In fact, our results are relatively insensitive to fsf. Although 
the wind dynamics do change during the “active” phase of the 
wind, most of the expansion occurs after active star formation 
has ended as long as fsf <C H{z)~^ ■ If this condition is satis¬ 
fied, the resulting bubble sizes do not depend strongly on fsf 
because bubbles spend most of their time in the late “remnant” 
stages. However, if fsf ~ //(z)“*, a significant amount of energy 
is lost to radiative cooling. In this case the total volume filled 
by metal-enriched winds can decrease by ~ 25%. 

With this prescription, we implicitly assume that a progenitor 
does not form stars and release mechanical energy until it joins 
the “trunk” of the merger tree. In reality, as long as the progen¬ 
itor mass exceeds Mniin(z), it would have formed stars before 
merging and thus be surrounded by its own wind bubble. The 
two bubbles will presumably merge along with the galaxies in 
some complicated manner. Our prescription is reasonable so 
long as the total wind energy is approximately conserved dur¬ 
ing the merging process. Note also that by choosing to follow 
only one of the progenitors at each branching, we do not nec¬ 
essarily track the largest single progenitor existing at each step. 
This error is not important so long as our Monte Carlo prescrip¬ 


tion includes a sufficiently large number of trials. 


2.2. Wind Luminosity 

Once we have fixed the rate of star formation /*Mg/fsf, the 
mechanical luminosity in the wind depends only on the ini¬ 
tial mass function (IMF) of the stars and the fraction of su¬ 
pernova energy available to power the wind. The IMF de¬ 
pendence can be characterized entirely by the total amount of 
stellar mass that must be formed in order to produce one su¬ 
pernova, wsN- For a Scalo (1998) IMF with lower and upper 
mass cutoffs M; = 0.1 Mq and M„ = 100 Mq, wsn = 126 Mq, 
while for a Salpeter (1955) IMF with the same mass cutoffs, 
WSN = 135 Mq. We use a Scalo IMF in the following because 
it better matches observations of the IMF in the local universe; 
however, the results are clearly insensitive to this choice. 

Recent simulations of the formation of metal-free stars indi¬ 
cate that these objects may be very different from nearby stars, 
with typical masses ~ 10^-10^ Mq (Abel, Bryan, & Norman 
2002; Bromm, Coppi, & Larson 2002). This mass scale is de¬ 
termined by the physics of H 2 cooling; we therefore assume 
that when molecular hydrogen is an active coolant, star forma¬ 
tion proceeds with a Very Massive Star (VMS) IMF. Unfortu¬ 
nately, these simulations do not yet have the dynamic range to 
follow the collapse to nuclear burning, so they are unable to 
predict the final IMF or the allowed mass range of metal-free 
stars. We therefore assume (somewhat arbitrarily) that the IMF 
is a Salpeter-like power law with an index (3 = 2.35 and mass 
cutoffs Ml = 100 Mq and M„ = 500 Mq. Interestingly, mas¬ 
sive, low-metallicity stars leave very different remnants from 
“normal” stars (Heger & Woosley 2002). For a stellar mass 
less than 130 Mq or greater than 260 Mq, the star collapses 
to a black hole without significant matter or energy ejection. 
However, stars with masses 130 Mq < M < 260 Mq are sub¬ 
ject to the e^e~ pair instability and explode without leaving a 
remnant black hole. Unlike normal supernovae, such explo¬ 
sions expel the highly-enriched core along with the envelope. 
This mechanism may thus have significant consequences for 
metal enrichment (see §3.1 below). For our VMS IMF, we find 
WSN = 462 Mq. 

We assume that each supernova releases an energy Fsn = 
lO^'ftsi erg into the interstellar medium of the host galaxy. 
For a Scalo IMF, we assume £51 = 1. For a VMS IMF, we 
assume £51 = 10 , an approximate logarithmic average of the 
mass-dependent explosion energy in the models of Heger & 
Woosley (2002). However, only a fraction /esc of this energy 
is available to power the wind; the rest is lost, primarily to ra¬ 
diative cooling, in the early stages of the outflow. Observations 
of the cool, dense ejecta of starbursts in the local universe re¬ 
quire /esc ^0.1 (Heckman et al. 2000), while observations of 
the warm and hot plasma in the wind imply /esc > 0.3 (Strick¬ 
land et al. 2000). However, at high redshifts the characteristic 
density is much higher, and radiative cooling could in princi¬ 
ple be more important. Simulations of high-redshift starbursts 
suggest that /esc ~ 0.25 (Mori, Ferrara, & Madau 2002). 

With these parameters, the mechanical luminosity driving the 
wind is 
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(3) 


2.3. Wind Model 
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We can now calculate the expansion of the wind into the 
IGM. The wind begins when the progenitor mass first exceeds 
Mmin(z), and energy is added to it during each merger event. 
We use the thin-shell approximation (Tegmark, Silk, & Evrard 
1993; Furlanetto & Loeb 2001), which is well-justified in a cos¬ 
mological context (Ikeuchi, Tomisaka, & Ostriker 1983). The 
expanding wind sweeps a fraction 1 - /m of the ambient IGM 
into a thin shell. The remaining gas leaks into the hot, rarefied 
interior whose pressure drives the shell expansion. Assuming 
spherical symmetry, we can then describe the outflow through 
the following system of equations: 

.. 4ttR^ G 

"'ur''’"'’’ 




P '■ 
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Here, R is the shell radius in physical units, Ms is the shell 
mass, and p is the pressure of the hot bubble interior. The shell 
expands through the host halo and into the IGM; the ambient 
baryon density, infall velocity field, and external pressure are 
denoted pg, Vinf, and pext, respectively. We assume that the host 
galaxy has photoionized the external medium, which is likely 
to be true during most of the wind’s expansion. The enclosed 
dark matter mass is and the host galaxy mass (which can 
be time dependent as the wind entrains material from the host) 
is Mgai. We use the universal halo profile of Navarro, Frenk, 
& White (1997) inside the accretion shock radius and the self¬ 
similar solution of Bertschinger (1985) outside of this radius, 
although the details of the profile are not critical to our results. 
Each halo is therefore embedded in a relatively small overdense 
infall region, which is in turn surrounded by IGM at the mean 
cosmic density. Both M/, and the halo profile evolve as the host 
galaxy grows through accretion. Finally, the rate of energy in¬ 
put L into the bubble interior is 


L — tsf)“t"Tdiss F 


comp 5 


(7) 


where Lcomp is the inverse Compton cooling luminosity, L„ is 
the wind luminosity defined in equation (3), 0 is the Heaviside 
step function, and 

L6Xss = -\fdM{R-VMf. ( 8 ) 


Here /j is the fraction of the energy dissipated by the accelera¬ 
tion of ambient material to the shell velocity that is retained by 
the outflow (rather than lost through cooling inside the shell). 
In the hot (T ~ 10^ K), rarefied bubble cavity, inverse Comp¬ 
ton cooling dominates at the high redshifts we consider, with a 
cooling time 


Comp 



(9) 


We therefore neglect metal line and free-free cooling inside the 
bubble."^ Material in the shell will cool rapidly through either 


inverse Compton cooling (when it is ionized) or line cooling 
(once recombination occurs). We therefore set fd = 0, although 
we also show some results for fd= I- We describe the expan¬ 
sion model in more detail in Furlanetto & Foeb (2001). 

As emphasized above, the structure of high-redshift winds 
is not well-constrained. They are likely to be more isotropic 
than those of local starbursts because their disks are not as well- 
defined. FBGs appear to host isotropic outflows, because wind 
absorption lines are seen in all FBGs that have been observed 
(Shapley et al. 2002). Moreover, the outflow will likely be¬ 
come more spherical as it plows into the IGM. We therefore 
assume that all winds expand spherically. A second key ques¬ 
tion is whether star formation occurs before or after the halo gas 
collapses into a dense medium. In the former case, the wind ex¬ 
pands into a spherical, relatively low density medium (see, e.g., 
Madau, Ferrara, & Rees 2001). In this “no disk formation” 
model, we assume that supernovae inject energy into the center 
of the halo and that all of the halo gas is swept up by the shell. 
Simulations show that a substantial fraction (^ 0.5) of the gas 
in the halo may remain inside the protogalaxy (Mori, Ferrara, & 
Madau 2002); however, this choice makes only a small differ¬ 
ence to our results because most of the expansion occurs after 
the swept-up IGM mass exceeds the halo mass (see §4.3). 

In most of our work, we assume that the halo gas collapses to 
a dense medium before the starburst begins, leaving the rest of 
the halo empty of baryons. We assume that the dense medium 
has a scale length 


^d ~ / —^vir; 


( 10 ) 


where R^„ is the virial radius of the halo as defined in equa¬ 
tion (24) of Barkana & Foeb (2001) and A = 0.05 is the spin 
parameter (Mo, Mao, & White 1998). With a slight abuse of 
terminology, we refer to this case as the “disk formation” sce¬ 
nario. Note that we do not require well-defined disks to form, 
only that the gas has collapsed to a dense central object. We 
assume that the wind entrains a fraction /s„ of the protogalac- 
tic material. Motivated by observations of local disk starbursts, 
we set /sw = 2/* in most of what follows (Martin 1999).^ The 
initial velocity of the outflow Vw is then fixed (c.f., Aguirre et 
al. 2001c; Springel & Hernquist 2002): 


Vw = 310 


/ /esc 2 126 Me \ 

\0.25 /sw//* WSN ) 


1/2 

km s“^ 


( 11 ) 


Note that Vw is independent of halo mass, implying that winds 
will most easily escape the gravitational attraction of smaller 
halos. The wind velocities in nearby starbursts are inferred to 
be independent of the host mass (Martin 1999). 

Our prescription for the initial conditions depends on 
whether disk formation is allowed. If it is not, we use the self¬ 
similar wind-fed solution of Ostriker & McKee (1988) to cal¬ 
culate the time for the wind to travel to Rd assuming a constant 
density within this region. (Although Rd has no physical mean¬ 
ing in this case, we use it to fix numbers.) The self-similar so¬ 
lution also yields the initial velocity, and we assume that the re¬ 
maining supernova energy goes into the bubble pressure. While 


^ Here we assume that the ejected metals are distributed in the same way as the shell and bubble gas, so that only a fraction /„ = 0.1 of the metals reside in the hot 
interior. Metal line cooling can then be safely ignored. If instead all the metals are confined to the hot cavity, metal line cooling cannot be ignored at z < 7 or when 
bubbles have cooled in the late stages of expansion (T < 10^ K); in other situations, inverse Compton cooling still dominates. At low redshifts, we note also that 
heating by the ultraviolet background radiation field is significant. 

^ Note, however, that the inferred star formation rates in some high-z galaxies (e.g., Pettini et al. 2001) are much larger than those of low-z starbursts. The calibration 
of Martin (1999) may therefore not apply to these systems. Fortunately, our results are relatively insensitive to this choice (see Figure 11). 
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this prescription makes several simplifying assumptions, the re¬ 
sults are insensitive to the precise initial conditions in this case. 
If disk formation is allowed, we assume that the wind travels 
with a constant velocity v„ within the dense collapsed region. 
The initial time is determined by the time taken for the wind to 
reach Rd- For numerical stability, we convert half of the kinetic 
energy to the thermal energy of the bubble interior before be¬ 
ginning the integration. (The final result is not sensitive to this 
conversion.) 

When gravity, density variations, external pressure, radia¬ 
tive cooling, and the cosmological expansion are ignored, our 
model for a single explosion recovers the usual Sedov-Taylor 
solution if fd= \ or the pressure-driven snowplow solution if 
fd = 0 (Ostriker & McKee 1988). With these effects included, 
we find instead R oc ^5 fj^ed halo mass, a somewhat steeper 
energy dependence than expected in self-similar solutions. This 
is principally due to the deceleration of the Hubble flow: an ex¬ 
plosion with smaller energy takes longer to escape the host halo 
and therefore must expend more energy accelerating the swept- 
up IGM. Our model matches well with previous analytic work 
on outflows (e.g., Barkana & Loeb 2001). However, it predicts 
wind radii approximately twice those found by Madau, Ferrara, 
& Rees (2001), even when we include only the effects that are 
part of their model. We do not consider this a cause for concern, 
because the simulations of Mori, Ferrara, & Madau (2002) have 
a similar discrepancy and compare favorably with our results. 

Figures 1 and 2 show the physical sizes R and expansion ve¬ 
locities Vexp of supernova wind bubbles at z = 20,12,8, and 3 
as a function of host halo mass. We assume the following pa¬ 
rameters: /* = 0.1, /esc = 0.25, wsN = 126 Mq (corresponding 
to a Scalo IMF), atomic cooling, disk formation, and fd = 0. 
The lower mass limit in each panel is simply M^in, and the up¬ 
per limit is the approximate mass M^ax above which the winds 
cannot escape the host halo. At a given mass, the dispersion 
in radius and velocity is caused by the randomly chosen star 
formation histories. For halos near Mmin, R varies significantly 
because such halos have undergone only a few star formation 
episodes. Thus, the time elapsed since the wind began fluctu¬ 
ates strongly. As M/, increases, the number of mergers increases 
and the dispersion decreases. Recall that the wind begins when 
the progenitor crosses the cooling threshold, and each subse¬ 
quent merger adds energy to it. Because massive halos have 
undergone many such random merger events, the fluctuations 
in merger times average out and the scatter in R is fairly small. 
The dispersion increases again nearMn,ax, where the star forma¬ 
tion history, detailed wind model, and initial conditions all have 
a substantial effect on whether the wind can escape beyond Ryir- 
Fortunately, halos near Mn,ax are sufficiently rare that they have 
negligible effects on our final results. 

R increases rapidly as redshift decreases because of the in¬ 
creasing cosmic time available to the outflow. The dispersion 
also increases with redshift, because small variations in the for¬ 
mation time make larger differences when the total lifetime of 
the outflow is small. The weak dependence of radius on mass 
is primarily due to the gravitational binding of the host. Be¬ 
cause the energy available to the starburst scales with while 
the gravitational potential scales with M^, we expect R{Mf) 
to flatten as mass increases until eventually the winds can no 
longer escape outside the virial radii of the hosts (shown by 
the solid curves). For most of the redshift range of interest 
and for reasonable choices of parameters, the threshold halo 
mass is M^ax ~ 10^®-3 x 10^' Mq, increasing as z decreases be¬ 


cause halos become less concentrated. Well below this thresh¬ 
old mass, we find R oc as would have been expected in 

a naive Sedov-Taylor solution (though the agreement is coin¬ 
cidental and the scaling only approximate). The top left panel 
of Figure 1 shows that LBGs can directly affect the surround¬ 
ing IGM to (physical) distances ~ 100 kpc, provided that they 
reside in halos with M/, ~ lO'^ Mq (see also Aguirre et al. 
2001b). Interestingly, Adelberger et al. (2002) find that regions 
within ~ 125 kpc of LBGs have unexpectedly small amounts 
of H I. One possible explanation is that winds have cleared this 
region of H I; if so, these observations suggest that our model 
provides a reasonable description of the extent of outflows. 


O 

D- 





Fig. 1.— Physical radius R of wind bubbles at various redshifts, as a func¬ 
tion of halo mass. The solid line in each panel shows the virial radius i?vir of 
the host halo. Note that both the mass and radius scales change between the 
different panels. All panels assume our standard set of parameters. 


Figure 2 shows that all of our wind bubbles have characteris¬ 
tic expansion velocities Vexp ~ 50 km s“', with a substantial dis¬ 
persion. The velocity increases slowly with mass, particularly 
at low redshifts, until M/, ~ Mmax- Our expansion velocities are 
much smaller than those inferred from observed starburst winds 
(Heckman et al. 2000; Pettini et al. 2002). The difference is, 
however, easily understood. While existing observations probe 
the early stages of starburst winds, our models primarily sample 
the late time behavior. At early times, the material ejected from 
the host galaxy dominates the dynamics, but gravity and col¬ 
lisions with the ambient IGM inevitably decelerate the winds 
later on. We can estimate Vexp by assuming that the host expe¬ 
riences one major merger per Hubble time and that all of the 
energy released goes into accelerating the swept-up medium. 
The average expansion velocity is then 


fexp 




126 Mq 
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10 

1 -i-z / 


Mh 

10^ Mq 
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( 12 ) 


This simple estimate neglects cooling, gravity, and the Hubble 
flow energy contained in the IGM, but it illustrates the approx¬ 
imate velocities we find in our model. Recalling R oc for 
Mh <C Mmax, equation (12) also predicts a weak mass depen- 
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dence Vexp oc , near what we observe in the complete model 
at small redshifts. Note that the Mach number J\4 = Vexp/Cj ~ 4, 
where Cj is the sound speed of the photo-ionized IGM. The thin 
shell model assumes ^ 1, so it is only marginally valid in 
this regime. The detailed expansion may therefore differ some¬ 
what from our model. 
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Fig. 2.— Expansion velocity Vexp of wind bubbles at various redshifts, as a 
function of halo mass. Note that the mass scale changes between the different 
panels. All panels assume our standard set of parameters. 

3. METAL LINES AT HIGH-REDSHIFT 

We now calculate the distribution of metals within the wind 
bubbles, including the most abundant elements and ionization 
states as well as the geometry of the metal-polluted regions. 

3.1. Metal Yields 

The relevant elements are easy to determine: they are those 
atoms most abundantly produced by Type II supernovae (for a 
Scalo IMF) or by pair instability supernovae (for a VMS IMF). 
We summarize the average yields per supernova for these el¬ 
ements in Table 1. For Type II supernovae, we use the yields 
of Woosley & Weaver (1995) averaged over a Scalo IMF with 
Z=1O-4Z0. These yields depend on the energy of the super¬ 
nova, which introduce uncertainties of ~ 25%. We have con¬ 
servatively taken the yields from the lower energy model. The 
metallicity of the progenitor star introduces a similar uncer¬ 
tainty, but we do not include chemical evolution within galaxies 
because other simplifications in our model create uncertainties 
at least as large. For the VMS scenario, we average the yields 
of Heger & Woosley (2002) over our power-law IMF. 

Table 1 also lists the properties of several interesting transi¬ 
tions, including the wavelength Am and the oscillator strength 
/osc. We will focus on these lines in the following. 

3.2. Wind Structure 

Next we must determine how the ions are distributed 
throughout the wind. Note that we assume purely spherical 
expansion throughout. In reality, each halo is surrounded by 
the “cosmic web” of filaments and voids. Dense filaments will 
slow the shell and break spherical symmetry. Fortunately, they 
fill only a small fraction of space and will not strongly affect the 


total volume of the bubble (e.g., Theuns et al. 2002). A more 
serious problem is fragmentation (see below). 

Observations of both local starburst winds (Heckman et al. 
2000) and LEG outflows (Pettini et al. 2002) show very com¬ 
plex ionization structure, with absorption from both low ion¬ 
ization states typical of H I regions (such as C II, O I, Si II, 
and Fe II) and high-ionization states (such as C IV and Si IV). 
Within the wind region, the hot gas is inferred to have a very 
high covering fraction Ca (defined as the fraction of the cross- 
sectional area covered by the gas), while the low-ionization 
gas has Ca ~ 0.4-0.8 depending on the source (Heckman et 
al. 2000; Shapley et al. 2002). 

If mixing between the ejected gas and the swept-up IGM is 
efficient, the thin-shell model requires that a fraction 1 -/m of 
the metals reside in a cool, dense shell with T ~ 10“* K, with 
the remainder in the hot, rarefied interior. Although there is 
good evidence for a shell structure in many observations of 
starburst winds (Martin 1998; Cecil et al. 2001), other obser¬ 
vations suggest that the cool gas resides in clouds entrained by 
the wind (Heckman et al. 2000). Even though there is often 
little evidence for a cool shell in starbursts of the latter type, 
the related observations exclusively study young starbursts in 
which the wind mass is still dominated by gas ejected from 
the host galaxy. Once the swept-up IGM gas dominates the 
wind, a shell will inevitably form (thought it may not be spher¬ 
ical). It is not clear whether the metals will fully mix with the 
swept-up material that constitutes most of the shell mass. How¬ 
ever, in our model the shell itself travels with Vexp ~ 50 km s“' 
while winds are launched with v„ ~ 300 km s“' (see equa¬ 
tion [11]). For M ^ Mmax, entrainment of the ambient IGM 
causes most of the deceleration; therefore metals in the wind 
will not decelerate significantly until they impact the shell at 
f ~ 5 X 10®(/?/kpc) yr, where mixing through turbulence or dif¬ 
fusion can occur. 

However, the shell may eventually fragment into clouds due 
to cooling, Rayleigh-Taylor instabilities, turbulence, or inho¬ 
mogeneities in the ambient medium. Fragmentation is partic¬ 
ularly likely given the relatively modest Mach number of the 
flow (see §2.3). At z = 3, an analysis similar to our own shows 
that galactic winds must fragment in order not to overproduce 
Lyct absorption lines (Theuns, Mo, & Schaye 2001). If ei¬ 
ther the shell has fragmented or the metals are distributed in 
cold clouds throughout the bubble, only a fraction Ca < I of 
all lines of sight through the wind bubble will produce metal 
absorption lines. On the other hand, each line of sight that 
does pass through a cloud will have on average a larger column 
density of metals. For unsaturated lines, the equivalent width 
W (X Nx (x C^', where Nx is the column density of the relevant 
ion. In this regime (which is satisfied for most of the lines of 
interest; see §4.1), we then have the transformation 

dN(>W) ^dNOCAW) 
dz dz 

where dN(> W)jdz is the number of absorption systems inter¬ 
sected per unit redshift with equivalent width greater than W 
along a random line of sight (see §4.3). For saturated lines, 
the equivalent width transformation factor will be smaller. Be¬ 
cause our results can be scaled relatively easily with Ca, for 
concreteness we assume that Ca = I- We further discuss the 
consequences of fragmentation in §5. 

The wind structure will also determine the width of the metal 
lines. In the shell (with T ^ 10“* K), the thermal line widths are 
bi<a ~ 3.1 {12!Ax)'^/^ km s“\ where b is the Doppler parameter 
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and Ax is the atomic mass of the ion x. In this case, though 
the lines may be very deep, they will be extremely narrow and 
hence difficult to detect except with extremely high-resolution 
infrared spectroscopy. However, low and moderate resolution 
spectroscopy of both nearby starbursts and LEG outflows show 
line widths comparable to or greater than the inferred mean ex¬ 
pansion velocity of the wind (Heckman et al. 2000; Pettini et al. 
2002). This can be easily understood if cloud fragments, each 
with an intrinsic line width bth, are distributed throughout the 
wind with a velocity dispersion comparable to Vexp, or if clouds 
are continuously produced and accelerated by the wind. We 
assume in the following that b = Vexp- Note that the expansion 
velocities we find (b ~ 50 km s“'; see §2.3) are much smaller 
than those typically observed in young starburst winds. If star 
formation is ongoing, and if the fast wind near the host contains 
a substantial fraction of the metals, we will underestimate the 
true line width. However, as long as the lines are not saturated, 
the equivalent width W is independent of our choice for b. Even 
if the lines are saturated, we cannot overestimate W so long as 
we underestimate the velocity dispersion. We emphasize again 
that the observations from which we calibrate b = Vexp use the 
starburst galaxy itself as the background light source; thus, they 
only constrain the active phases of the starburst. The old winds 
that we study through background quasar spectra may have very 
different velocity dispersions. 


3.3. Ionization States 


Given the most abundant elements, we now determine the 
relevant ionization states. Three ionizing processes must be in¬ 
cluded; collisional ionization (especially in the hot interior), 
photoionization from the host galaxy, and photoionization by 
the extragalactic background light. In general, if a fraction Cm 
of a given element is in a specified ionization state, the covering 
fraction of the outflow is unaffected but Nx (x Cm- Thus in the 
unsaturated line limit 


dN(> W) dN(> C~‘W) 
dz dz 


(14) 


Cm < 1 can also occur if metals are retained by the host galaxy. 
However, simulations suggest that the wind carries a very large 
fraction of the metals produced during a starburst (Mac Low & 
Ferrara 1999). 

In the thin shell model, we can estimate the ionization state 
of the metals in each phase. We first consider metals in the 
shell. If the hydrogen in the shell is neutral, metals will reside 
in states typical of H I regions (C II, O I, Si II, and Fe II for 
the elements of interest), because photons able to further ionize 
these elements will be absorbed by H I. The ionization poten¬ 
tial of OI is 13.62 eV, and it remains locked in charge exchange 
equilibrium with H I (see Oh 2002 for a discussion of the util¬ 
ity of this mechanism). The other elements have first ionization 
potentials below that of hydrogen. We now show that the extra- 
galactic background light most likely suffices to maintain these 
first ionized states. The total column density of a metal x along 
a typical line of sight through (each side of) the wind is 


Yx fib f,Mh 

WSN ^0 ATTAxMpR'^ ' 


(15) 


We can estimate the background flux /(t'ion) = Tion, where r'ion is 
the ionization threshold of the relevant species, needed to keep 
a metal ionized by requiring that the ionization rate exceed the 
recombination rate within the shell: Jion ^ OixngNx.iot, where 


is the recombination coefficient and n^ is the electron density in 
the shell. We find 

126 Mq Alfi 
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in units of erg s"' cm-^ Hz"' sr"* and where hb is the mean 
baryon density of the IGM. The electron density is difficult 
to constrain. The standard jump conditions for an isothermal 
shock give a shell baryon density Wsh ~ Ai^Abb, where A is the 
local overdensity through which the shell travels. However, the 
shell initially plows through highly overdense material within 
the halo (A > 200) at a very high velocity. Later, the shell 
thickness can only increase at the sound speed and it may not 
have had sufficient time to expand to its equilibrium structure. 

We now estimate the extragalactic background over the range 
~ 7-13 eV, which includes photons able to ionize C, Si, and Fe. 
We assume that the emissivity of a typical galaxy obeys e oc i'~°‘ 
and is proportional to the mass of the galaxy. The volume emis¬ 
sivity of photons is then 


j{v-,z) = 2'KafiNu 


af*n% 1 



dz 

^^coll 

\^h) 

dt 

dz 


(17) 


where Fcoii is the fraction of collapsed baryons in galaxies and 
where we have normalized the emission spectrum to the num¬ 
ber of photons per baryon in stars emitted above the ionization 
threshold of hydrogen i^h, Mon ~ 4000 (Barkana & Loeb 2001). 
The extragalactic background is (Haiman, Rees, & Loeb 1997) 

JAz)=^ [ j(izz,z')e~'"^"^'^d£(z'), (18) 


Att 


where = v[{l+z')/(l+z)], t(z') is the IGM optical depth at 

and £ is the proper length. Below iz/j, the only sources of signif¬ 
icant opacity are the Lyman lines of hydrogen. For a given ob¬ 
served frequency iz at z, the universe is thus optically thin back 
to the point Zon at which this radiation would have blueshifted 
into the nearest Lyman line (denoted here as Lyx), 


1+z, 

Thus (in units of erg s"* 
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( 20 ) 


For carbon, the integral is typically (~ 10' -3,3 X 10-4,2 X 
10-3,10-3) for z = 20,16,12, and 8, with weak dependence on 
the spectral slope. The background for Fe and Si is a few times 
larger. Comparing equations (16) and (20), we see that so long 
as the shell is neutral (n^ <C Ush) or the shell is not extraordinar¬ 
ily overdense, the extragalactic background will easily ionize 
these elements to C II, Si II, and Fe II. 

By definition, the extragalactic background cannot ionize hy¬ 
drogen in the shell until after cosmological reionization. We 
now consider whether the host galaxy can ionize the shell. 
According to the starburst models of Leitherer et al. (1999), 
the number of ionizing photons produced per unit time by an 
instantaneous starburst with a Salpeter IMF and metallicity 
Z = O.O5Z0 (the lowest metallicity they calculated) a time t after 
the burst is 

^■^7 ^ ^ ^ \ g-l 

dt Vio’yry 


(21) 



provided that t>10^ yr. Here /*Mg is the stellar mass produced 
in the starburst. Of these, only a fraction x can escape the host 
galaxy. This escape fraction is poorly constrained at present 
but is likely < 10% (Leitherer et al. 1995; Wood & Loeb 2000; 
Dove, Shull, & Ferrara 2000; Steidel, Pettini, & Adelberger 
2001). 

As above, we can calculate the column depth of ionized hy¬ 
drogen Ahii by balancing the number of ionizations and recom¬ 
binations per second. The recombination rate again depends on 
the shell density. Ahu will be largest for the smallest we 
therefore assume pessimistically that Wsh ~ M.^nb- We com¬ 
pare Ahii to the total column in the shell Ash, assuming that the 
swept up matter dominates the shell mass. We then find, for an 
instantaneous starburst, that 
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We have confirmed the approximate validity of this simple 
estimate over the range of densities and spectra of interest 
with CLOUDY (Ferland 2001). Thus, once an appreciable 
time has elapsed after the starburst, the shells remain neu¬ 
tral. However, if active star formation persists, dN^/dt ~ 
M© yr) s“' (Leitherer et al. 1999). The shell could 
then in principle remain partially ionized for a longer period. If 
the shells fragment or if the metals otherwise reside in clouds, 
the densities will increase beyond Ai^nn, making ionization 
even more difficult. In this case, the ionization balance would 
also depend on the size of the clouds, parameterized in our 
model by Ca- 

Observations of nearby starburst winds and LEG outflows 
find strong absorption from low ionization states (Heckman et 
al. 2000; Pettini et al. 2002). In all of these systems, the star- 
burst is still active and hence dN^ jdt is large. In the old winds 
that provide most of the absorption lines, the arguments above 
show that low ionization states should be even more important. 
In the following we therefore focus on low ionization states, 
and we assume that Cm = 1. While this is obviously a naive 
approximation, our results can easily be rescaled by equation 
(14). For completeness, we also show some results for ions ap¬ 
propriate to a “warm” medium, C IV and Si IV. These ions may 
be appropriate if metals are advected to the shell but do not mix 
efficiently with the swept-up IGM, thus remaining at the inter¬ 
face between the cool shell and the hot cavity. Lines from these 
ions are prominent in LEG winds (Pettini et al. 2002; Shapley 
et al. 2002 ). 

If, however, the metals remain locked in the hot bubble cav¬ 
ity, our assumptions will break down. This medium has density 
~ fniiH and is easily ionized by even trace amounts of star for¬ 
mation in the host. More importantly, the gas in this phase has 
typical temperatures T ~ 10^ ^-10® ^ K and thus is collision- 
ally ionized to very high ionization states that lack rest-frame 
optical or ultraviolet absorption lines. Fortunately, any metals 
carried by the hot wind medium will eventually be advected to 
the shell, where cooling and mixing will occur (see §3.2). 


4. RESULTS 

With the wind model described in §2 and the metal charac¬ 
teristics outlined in §3, we can now calculate the consequences 
of metal pollution from galactic winds. We list the parameters 


of our model, along with their standard choices, in Table 2. Un¬ 
less otherwise specified, all curves shown assume these values. 
We also consider two halo cooling channels: atomic hydrogen 
cooling, which allows halos with Tvir > 10"^ K to form stars, 
and molecular hydrogen cooling, which reduces the threshold 
to Tvir > 400 K. In the former case, we assume a Scalo (1998) 
IMF applies, while in the latter we assume a VMS IMF. This 
choice in turn determines wsn, Es\, and Y^. For H 2 cooling, we 
set wsN = 462 M© and £51 = 10. 

We present the results of our model in the following sections. 
First we describe the observable characteristics of individual 
winds in §4.1. We then discuss the global evolution of the metal 
filling factor in §4.2. Finally, we present our predictions for ab¬ 
sorption line statistics in §4.3. 
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Fig. 3.— Central optical depth to of the O I A1302 feature as a function 
of halo mass for our standard parameter choices. Crosses, horizontal lines, 
triangles, and points refer to z = 20,16,12, and 8, respectively. 



4.1. Wind Characteristics 


Absorption lines are characterized by three basic parameters 
(of which two are independent): the line width b, the central 
optical depth tq, and the equivalent width W. In our model, we 
set b = Vexp (see Figure 2). Although we present our results ex¬ 
clusively in terms of W, it is also useful to examine the central 
optical depth tq of the absorption systems. In a real observa¬ 
tion, To helps to determine the ease with which each line can 
be identified. It also indicates where saturation becomes im¬ 
portant (and hence where equations [13] and [14] break down). 
The typical column density of an ion x through (each side of) 
the wind per velocity interval is 


= (23) 

s/TTb 

where u is the velocity relative to Vexp- The optical depth at line 
center is then 


TO = 


TreVoscA 

m N, ,tot 

MgC 


(24) 


f /osc Y, 

16 A„ A 

y50kms-‘^^ 

\0.05 0.5 M© 

Ax 1300 Ay 

1 b ) 


/* 126 M© Oi/Uo Mb A /20kpc\ 
^ wsN 0.05/0.3 10® M ©) R ) 




Furlanetto & Loeb 


9 


We show To as a function of halo mass for several redshifts in 
Figure 3. We have chosen the O I A1302 transition for this plot 
because it is the strongest low ionization absorber in our model. 
Because R ex while b is roughly constant in our model, 

To (X below Mniax- To increases with z because wind bub¬ 
bles are more compact at higher redshifts and hence the typical 
column density Nx is much larger. We see that line saturation 
is not important for the low mass-halos that contribute most of 
the filling factor at z < 12, but it cannot be neglected at higher 
redshifts. At a given z, the minimum to is determined by the 
smallest halos that can collapse to form stars (fixed in this case 
by the atomic cooling threshold). 
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Fig. 4.— Observed equivalent width W of the O I A1302 feature as a func¬ 
tion of halo mass M/, for our standard parameter choices. Crosses, horizontal 
lines, triangles, and points refer to z = 20,16,12, and 8 , respectively. 
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(Kauffmann et al. 2002; Dekel & Woo 2002). If so, this would 
steepen the Mh-W relation while preserving the one-to-one cor¬ 
respondence. 



Fig. 5.— Filling factor of polluted regions as a function of redshift. Solid 
curves assume H 2 cooling and a VMS IMF, while dotted curves assume atomic 
cooling and a Scalo IMF. Each set has /* = 0.5,0.1, and 0.01 from top to bot¬ 
tom; other parameters take their standard values. The dashed curve assumes 
H 2 cooling with /* = 0.1 and £’51 = 1 . 


4.2. Filling Factor 


To calculate the filling factor of the wind bubbles, we com¬ 
pute 


Qiz) = 




dM^V{M,z), 
dm 


(27) 


The observed equivalent width W is 

W = (l+z) j [l-e-^<^>] d\. 

If the line is not saturated, this expression may be written 
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We show W as a function of halo mass in Figure 4 (again for 
O I A1302). Typical equivalent widths increase from ~ 0.5- 
10 A between z = 8 and z = 20. W increases with redshift both 
because of stretching by the cosmological redshift and because 
the wind bubbles are smaller at these times. When line satura- 
tion can be neglected, W oc bTo oc . As the lines become 
saturated, the relationship flattens. Interestingly, Figure 4 sug¬ 
gests that a line of a given equivalent width can be associated 
with a halo in a fairly narrow mass range, particularly at lower 
redshifts. This is essentially a consequence of the relatively 
small dispersion in wind bubble sizes for a fixed set of param¬ 
eters and our assumption that /* and /esc do not vary systemat¬ 
ically with halo mass. As a result, Nx is well-determined given 

2/3 

the halo mass. There are some recent indications that /* oc Mr 


where V(M,z) is the comoving volume of the wind produced 
by a halo of mass M at a redshift z and dn/dM is the halo mass 
function (Press & Schechter 1974; Sheth & Tormen 1999). 
Note that we do not include photoheating due to reionization 
in determining Mmin. We therefore likely overestimate Q{z) for 
z < 6, but we include these results for comparison with previous 
work. Because the star formation history determines the wind 
size, we perform the integral using a Monte Carlo algorithm. 
We chose the mass sampling by requiring that our results con¬ 
verge to < 10% between trials.® Because winds escape small 
halos most easily, Q{z) is dominated by winds from low-mass 
halos (in most of our calculations, winds with Mh < lOMnjin(z) 
account for ~ 75% of the total polluted volume); in practice, the 
integral need only extend to some maximum mass M^ax (see 
§2.2). The principal shortcoming of this model is that it ne¬ 
glects the overlap of wind bubbles; thus, we can have Q{z) > 1. 
If sources are randomly distributed, the true filling factor will be 
$ = 1 -exp[-2(z)]. However, the relationship between Q{z) and 
the true filling factor cannot be so easily quantified if sources 
are clustered (see §5). 

Figures 5 and 6 show Q(z) in our model. In Figure 5, we 
compare Q(z) for H 2 cooling (solid curves) and atomic cooling 
(dotted curves). Within each set, the three curves correspond to 
/* =0.01,0.1, and 0.5, from bottom to top. All other parameters 
have their standard values. We see that, in the atomic cooling 


® For ftfese ^5x10^, stochastic fluctuations are iarge and the convergence is somewhat worse; see Figure 5 for an exampie. 
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scenario, Q(z = 12) « 1% and Q(z = 6 ) « 10% (for /* = 0.1). 
The filling factor increases rapidly with cosmic time because 
of the increasing fraction of baryons in collapsed objects and 
because the winds have considerably more time to expand at 
low redshifts, a point emphasized by Aguirre et al. (2001c). By 
allowing less massive halos to form stars and by injecting an 
order of magnitude more energy per supernova, the H 2 cool¬ 
ing scenario leads to a much larger filling factor: with /* = 0 . 1 , 
this model has Q(z = 12) « 1. The dashed line (nearly obscured 
by the lowest solid line) shows the H 2 cooling scenario with 
/* = 0.1 and £51 = 1 ; thus, even if H 2 cooling creates “normal” 
massive stars the enrichment process is still much more effi¬ 
cient than with pure atomic cooling. At very high redshifts, 
the minimum halo mass in the atomic cooling scenario is above 
the nonlinear mass scale in the Sheth & Tormen (1999) mass 
function. In this regime, the halo abundance is exponentially 
suppressed; therefore, Q(z) for H 2 cooling greatly exceeds that 
for atomic cooling. The gap closes (though it is still substan¬ 
tial) as Mniin falls below the nonlinear mass scale at lower red¬ 
shifts. Our results for Q(z) in the atomic cooling scenario are 
very similar to those of Scannapieco, Ferrara, & Madau (2002), 
even though they use a more sophisticated algorithm to deter¬ 
mine star formation efficiencies and halo clustering. We find 
significant differences, however, between their “Population III” 
scenario and our H 2 cooling scenario. The disparity is likely 
caused by their requirement that H 2 cooling occurs before disk 
formation which imposes a much more stringent mass thresh¬ 
old than ours, and because we assume that VMS supernovae are 
more energetic. 



Fig. 6 .— Filling factor of enriched regions as a function of redshift normal¬ 
ized to /*/esc, 6 (z)/(/* /esc). The three solid curves show /esc = 1,0.25, and 
0.1 with /* = 0.1. The dashed and dotted curves assume/* = 0.5 and /* = 0.01, 
respectively, with /esc =0.25. All other parameters take their standard values. 


Figure 6 plots Q(z) /(/*/esc) for the atomic cooling scenario. 
We have thus normalized Q(z) by the factor controlling the en¬ 
ergy output in the outflows. This normalized filling factor is 
constant (to within a factor ~ 2 ) at a fixed redshift as long as 
/*/esc ^ 5 X 10“^ (c.f., Scannapieco, Ferrara, & Madau 2002). 
As discussed in §2.2, our model yields R oc £’•' 35 ^ fixed mass 

Ml, <C Mmax- Thus, 0(z)/(/*/esc) should be nearly constant so 
long as the halos that contribute most to the filling factor (i.e.. 


the first decade in mass above M^nin) are well below this mass 
threshold for efficient escape. The proportionality is broken for 
/♦ /esc ^ 5 X 10“^ because the number of halos hosting substan¬ 
tial winds begins to depend on /*/esc. 


4.3. Line Statistics 


Unfortunately, the volume filling factor Q(z) is difficult to 
observe, requiring sophisticated pixel-by-pixel analyses (e.g., 
Aguirre, Schaye, & Theuns 2002). It is therefore more useful 
to consider the statistics of absorption lines. The number of sys¬ 
tems intersected per redshift interval with an observed equiva¬ 
lent width W greater than some threshold Wo along a random 
line of sight is 


dN(> Wo) 
dz 


dr 

dz 



dM^A(M,z)e(W-Wo), 

dm 


(28) 


where A(M,z) is the comoving cross-sectional area of the wind, 
dr/dz is the comoving distance per unit redshift, and 0(W-Wo) 
is unity if the equivalent width of a typical line of sight through 
the wind is greater than Wo and zero otherwise. Again, our 
model neglects halo clustering (see §5 for a discussion of its ef¬ 
fects). We could in principle choose random impact parameters 
for each bubble; however, given the approximations we have 
already made, we consider this extra layer of complexity to be 
unnecessary. Note that equation (28) assumes a single absorp¬ 
tion line per halo. In the ideal case, each intersected wind bub¬ 
ble would cause 2 lines separated by ~ 2vexp. However, if the 
velocity resolution of the observation is lower than this value, 
the two lines would appear as a single feature, albeit with an 
equivalent width up to twice that which we assume. The two 
lines will also merge if b > Vexp- For clarity, we therefore as¬ 
sume that only one line is produced per halo. 

Another important quantity is a measure of the total flux 
per redshift interval absorbed by a single transition (neglecting 
overlap of the lines). 


Wtot = 


/wA 

'dN(> wy 

/ dW 

dz 


dW. 


(29) 


If Wtot Am, lines will be isolated and easily separated. If 
Wiot ~ Am, lines blend together and absorb a large fraction of 
the flux from the background source. 

Figure 7 shows the evolution of dN(> W)/dz with redshift. 
Here we show results for O I A1302, the strongest line in the 
atomic cooling scenario. The curves show line statistics at 
z = 20 (dotted line), z = 16 (short-dashed line), z = 12 (long- 
dashed line), and z = 8 (solid line). Each assumes our stan¬ 
dard parameter set. Clearly, the number of lines intersected de¬ 
pends strongly on redshift, particularly for weak lines (increas¬ 
ing from ~ 0.25 to ~ 40 with W > 0.5 A per redshift interval 
from z = 20 to z = 8 ). Two factors contribute to this increase: the 
winds are more compact at higher redshifts (see Figure 1) and 
the fraction of baryons in star-forming galaxies is much smaller 
at higher redshifts. The shape of the curves are easy to under¬ 
stand. As shown in Figure 4, a given W roughly corresponds 
to halos of a particular mass at fixed z. Thus, at each redshift, 
dN(> W)/dz is a measure of the (cumulative) halo mass func¬ 
tion weighted by R^. The “saturation” at small W is simply a 
result of the minimum halo mass for star formation Fig¬ 
ure 7 also shows how Wtot evolves with redshift. Recalling that 
OI is the strongest line we consider, line blanketing should not 
be significant for /* = 0.1 in the atomic cooling scenario. 
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Fig. 7. — Number of intersections per redshift interval above a given equiv¬ 
alent width threshold for O I A1302. From top to bottom, curves refer to 
z = 8,12,16, and 20. All parameters take their standard values. We also list 
Wtot (see equation [29]) for each line. 

Figure 8 shows the O I statistics for the atomic cooling sce¬ 
nario at z = 8 . From top to bottom, the three solid curves as¬ 
sume /h, = 0.5,0.1, and 0.01, with /esc = 0.25. The two dashed 
curves assume /esc = 1 and 0.1 (top and bottom at the left, re¬ 
spectively), with /* = 0.1. All other parameters are standard. 
Recall that Q(z) oc /*/esc; clearly, the dependence of the line 
statistics on these parameters is not degenerate. Essentially, 
absorption lines give information about both the extent of the 
enriched IGM and the level of enrichment. Recall that /esc af¬ 
fects only the amount of energy powering the wind, while /* 
affects both the energy input into the wind and the total mass of 
metals produced. Thus, an increase in /esc increases the num¬ 
ber of weak lines but decreases the number of strong lines, be¬ 
cause the metals become more dilute. An enhanced star for¬ 
mation efficiency increases the number of lines of all strengths. 
The total absorption Wtot strikingly reveals the difference. As 
shown in Figure 8 , Wtot increases rapidly with /*. However, it 
varies only slowly with /esc, increasing from 41 A to 57 A be¬ 
tween /esc = 0.1 and /esc = 1. Note that if /* is very high at 
this redshift, line blanketing from the many possible transitions 
becomes important. 

The greatly increased filling factor of polluted regions for 
H 2 cooling compared to atomic cooling (see Figure 5) suggests 
that these two scenarios may be distinguished through obser¬ 
vations of line statistics. Figure 9 shows that this is indeed 
the case. We show the statistics for four different transitions 
(O I A1302; long-dashed lines; C II AI334, solid lines; Fe II 
A2344,^ dotted lines; Si II A1304, short-dashed lines) at z = 20 
and z = 12. The upper set of curves in each panel assumes H 2 
cooling and a VMS IMF, while the lower set assumes atomic 
cooling and a Scalo IMF. All other parameters are standard. 
We also list Wot for the O I A1302 transition in each case. The 
difference between dN/dz for the two models (a factor ~ 1100 
at z = 20) is even more dramatic than the difference in Q{z) 
(a factor ~ 500 at z = 20). This occurs because of the rapidly 


increasing mass function at small masses, which allows many 
more individual halos to produce winds if H 2 cooling is permit¬ 
ted; the line statistics are more sensitive to the number of star¬ 
forming halos than is Q{z) because they depend most strongly 
on the total number of individual winds that can in principle be 
intersected. However, because the additional halos are small, 
they generally produce weak lines. Therefore the enhancement 
of weak lines is much larger than the enhancement of strong 
lines (this being particularly evident at z = 12). Interestingly, 
Wot is very large for all redshifts in the H 2 cooling scenario. 
Thus a large fraction of the background source flux could be 
absorbed by overlapping metal lines in this scenario, extending 
the apparent Gunn-Peterson absorption to much higher wave¬ 
lengths than otherwise expected. 



Fig. 8 .— Number of intersections per redshift interval above a given equiv¬ 
alent width threshold at z = 8 for the O I A1302 transition. The solid lines 
assume /* = 0.5,0.1,0.01, from top to bottom, and /esc = 0.25. The dashed 
lines assume /esc = 1 and 0.1, from top to bottom at left, and /* = 0.1. All 
other parameters are standard. We also list Wtot (see equation [29]) for the 
solid curves. Wtot varies only slowly with /esc (see text). 

We also see the effects of a VMS IMF in Figure 9. Although 
WSN varies by less than a factor ~ 4 between the two models 
(and is actually larger for a VMS IMF), the explosion energy 
for the VMS IMF is an order of magnitude larger than for a 
Scalo IMF. (Note, however, that even if £51 = 1 for VMS stars, 
the differences between the two scenarios are large; see Figure 
5). The different yields of the two IMFs also cause differences 
in the relative line strengths. For a Scalo IMF, the strongest fea¬ 
tures are O I and C II, simply because these are the most abun¬ 
dant elements produced in Type II supernovae. Pair instability 
supernovae, on the other hand, completely disrupt the star and 
therefore expel heavy elements much more efficiently (Heger & 
Woosley 2002; see Table 1). As a result, the strongest features 
in VMS winds are O I and Fe II, with C II being the weakest 
feature. 

We argued in §3.3 that the low ionization states of metals are 
likely to be most important, but it is also possible that many 
metals will be in relatively high ionization states. Figure 10 
shows dN(> W)ldz at z = 8 . The thick curves show low ion- 


’ We have selected this line from the several available Fe II transitions because it has an “average” oscillator strength. The strongest line, Fe II A2382, will have an 
equivalent width ~ 3 times greater than that shown. 
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ization states (as in Figure 9) while the thin curves show high 
ionization states (C IV A1548, solid line; Si IV A1394, short- 
dashed line). In each case we assume Cm = 1- We see that high 
ionization states actually provide a slightly stronger signal than 
the low ionization states; our previous predictions are there¬ 
fore conservative in this sense. However, if the metals are even 
more highly ionized, no observable transitions will remain with 
Am ^ Alyo ■ 



Fig. 9. — Number of intersections per redshift interval above a given equiv¬ 
alent width threshold. In each panel, the bottom set of curves assumes atomic 
cooling while the top set allows H 2 cooling. All other parameters take their 
standard values. We also list Wiot (see equation [29]) for the OI A1302 transi¬ 
tion in each scenario. 

Finally, Figure 11 shows the effect of varying other param¬ 
eters of the outflow model. All the curves assume our stan¬ 
dard parameter choices with the following exceptions: fd = l 
(dotted line), /sw = 1 (dot-dashed line), and no disk formation 
(dashed line). We see that the uncertainty due to model details 
at z = 8 is < 30%. The uncertainties in Q{z) are somewhat larger 
(~ 50%), because the integrand of equation (27) is proportional 
to R^, while dN/dz is proportional to the cross-sectional area 
R^. Thus our predictions for line statistics are fairly insensitive 
to the detailed wind parameters. Note, however, that the uncer¬ 
tainties increase with redshift, reaching a factor ~ 2 by z = 20 
for atomic cooling. At high redshifts, the winds have spent a 
larger fraction of their expansion within the halo and are hence 
more sensitive to its structure. In contrast, at low redshifts most 
of the expansion takes place in the IGM and is independent 
of assumptions about the halo structure. Model assumptions 
are much less important in the H 2 cooling scenario because the 
characteristic halo mass (and hence physical scale over which 
the halo structure matters) is much smaller. 

We have examined the absorption by metals in high-z star- 
burst winds seen against a background light source. We allow 
each star-forming halo to host a wind, with a time-dependent 
luminosity determined by its star formation history. We use 
the thin-shell approximation to trace the size of each wind bub¬ 
ble. Our model predicts that LEG winds can reach a proper 
distance ~ 100 kpc from the host. It is encouraging that this re¬ 


sult roughly matches the observationally inferred distance from 
LBGs within which HI is suppressed (Adelberger et al. 2002).* 
Our results for the filling factor Q(z) are consistent with previ¬ 
ous work (e.g., Scannapieco, Ferrara, & Madau 2002), but we 
have emphasized the large uncertainties in the wind modeling 
and the correspondingly large uncertainty in Q(z). We argue 
that the statistics of metal absorption lines are a useful way to 
constrain these parameters. Not only are these statistics directly 
observable, but they are also less dependent on the wind model¬ 
ing than Q(z) is. Our predictions are therefore relatively robust, 
although detailed modeling of the wind structure will be essen¬ 
tial for a complete interpretation of the line statistics. 



Fig. 10. — Number of intersections per redshift interval above a given 
equivalent width threshold at z = 8, assuming our standard parameter choices. 
Thick solid, long-dashed, short-dashed, and dotted curves are for C II, O I, 
Si II, and Fe II, respectively. Thin solid and short-dashed curves are for C IV 
and Si IV, respectively. In each case we fix Cm = 1. 

5. DISCUSSION 

We estimate that ~ 40 OI A1302 lines with equivalent width 
W > 0.5 A should be visible per unit redshift at z = 8, assuming 
a star formation efficiency /* =0.1 and neglecting fragmenta¬ 
tion. This decreases to ~ 0.25 lines per unit redshift by z = 20, 
both because of the smaller fraction of baryons in collapsed ob¬ 
jects and because of the brief time that each wind has had to 
expand by that time. The normalization and shape of the line 
statistics will primarily constrain /* and /esc, the fraction of su¬ 
pernova energy that powers the wind. Because the total amount 
of absorption is proportional to the amount of metals ejected, it 
depends strongly on /* but is nearly independent of /esc. The 
ratio of strong to weak lines, on the other hand, measures /esc 
because that parameter controls the volume over which a fixed 
metal mass can expand (see Figure 8). 

If we allow H 2 cooling to occur, then we expect ~ 250 O I 
A1302 lines per unit redshift at z = 20, because this cooling 
mechanism lowers the threshold halo mass for star formation 
and increases the explosion energy per supernova. In fact, the 
combined absorption from the several allowed transitions will 
eliminate a substantial fraction of the flux from the background 


^ Note that we do not include photoevaporation or infall suppression after reionization (see Barkana & Loeb 2001 and references therein) in our model star formation 
histories. However, LBGs are the most massive objects at z ^ 3, for which these processes can be ignored. 
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Fig. 11.— Number of intersections per redshift interval above a given 
equivalent width threshold at z = 8 for the OI A1302 transition. All curves as¬ 
sume our standard parameter set with the following exceptions: = 1 (dotted 
curve), no disk formation (dashed curve), and /sw = 1 (dot-dashed curve). 

source at wavelengths longer than the onset of the Gunn- 
Peterson trough. This would be straightforward to detect, but 
would also make isolating individual lines difficult. The dra¬ 
matic differences between this scenario and the atomic cooling 
scenario are large enough that the two should be easily separa¬ 
ble even given our uncertainties in /*. Metal absorption lines 
therefore offer a possible route to determine when and if H 2 
was destroyed throughout the IGM. This transition has impor¬ 
tant implications for the early growth of cosmic structure (e.g., 
Wyithe & Loeb 2002d; Mackey, Bromm, & Hernquist 2002). 
The relative strengths of the observed lines contain information 
about the IMF. In particular, Fe lines will be much more promi¬ 
nent if pair instability supernovae dominate (i.e., if the IMF is 
biased toward very massive stars). 

Metal absorption line studies can shed light on the chemi¬ 
cal evolution of the universe. Qian & Wasserburg (2001, 2002) 
have proposed an important model for the enrichment history 
of the universe based on abundance patterns in Galactic halo 
stars. They argue for a “prompt” enrichment phase to a thresh¬ 
old iron abundance [Fe/H] « -3 after which the abundance 
patterns change markedly; below this threshold, VMS super¬ 
novae dominate the enrichment, whereas above it “normal” su¬ 
pernovae dominate. The authors implicitly assume, based on 
ionic abundances in the Lya forest (Qian et al. 2002), that the 
prompt phase uniformly enriches all regions of the IGM to this 
threshold. Qian et al. (2002) find little evidence for evolution 
in the elemental abundances in the Lya forest over the range 
z ~ 4.6-0.09, indicating that the prompt enrichment occurred 
at z > 5. Qualitatively, we have shown that an early phase 
of halo formation through H 2 cooling forming VMS stars can 
pollute large regions of the IGM, lending credence to this as¬ 
sumption. Figure 5 shows that complete enrichment of the en¬ 
tire IGM by this epoch requires fairly efficient star formation 
if* ^ 10%) in small halos. In fact, observations of the Lya for¬ 
est metallicity do not require complete mixing of metals in the 
IGM, because these studies are currently only sensitive to ab¬ 
sorbers with Nhi ^ 10 ''^ ® cm“^, corresponding to moderately 


overdense regions. Metal enrichment is likely to be confined 
in islands that correlate with overdense regions where the Lya 
absorption is enhanced. If VMS regions are preferentially clus¬ 
tered near cosmic overdensities, it is plausible that those re¬ 
gions observable as Lya absorbers with HI column densities 
> 10 '^ ^ cm“^ are enriched while other regions remain pristine. 
In this case, lower star formation efficiencies would be allowed. 
We would also expect regions polluted by VMS supernovae and 
those subsequently polluted by normal supernovae to differ in 
their abundance patterns. An observational probe of the abun¬ 
dance patterns in the low-density IGM could therefore reveal 
the era over which the pollution occurred, although any such 
measurement would be extremely difficult. Note that the Qian 
& Wasserburg (2001) model assumes a one-to-one, monotoni- 
cally increasing relation between metallicity and cosmic time. 
If mixing is not complete at high redshifts, the transition from 
“very massive” to “normal” star formation may occur at differ¬ 
ent times in different regions of the universe. Thus, the non¬ 
uniformity of the enrichment process, which metal absorption 
line studies probe directly, determines the non-uniformity of the 
transition between different modes of star formation. 

We argued that a substantial fraction of the metals are likely 
to exist in low-ionization states (typically C II, O I, Si II, and 
Fe II). If, however, the metals exist primarily in higher ioniza¬ 
tion states, such as C IV and Si IV, then the absorption signal 
will actually be somewhat easier to observe because the rele¬ 
vant transitions are stronger (see Table 1). Moreover, each of 
these transitions is a doublet and hence relatively easy to iden¬ 
tify unambiguously. If the metals remain in the hot bubble gas 
with T ~ 10® K, they will reside in extremely high ionization 
states (such as O VI) whose resonant transitions are all blueward 
of Lya and therefore unobservable in the high-z IGM. 

The largest uncertainty in constraining /*, /esc, H 2 cooling, 
and the IMF will probably be the covering factor of the enriched 
material in the wind, Ca- Meaningful constraints on this fac¬ 
tor must await more detailed models of high-z winds, because 
inferences from observations of local starburst winds may not 
apply to sources in the early universe. For example, there are in¬ 
dications that outflows from LBGs are approximately spherical 
(Shapley et al. 2002), unlike the bipolar winds observed in local 
starbursts. Furthermore, existing observations of absorption by 
starburst winds and LBG outflows all use the host galaxy as the 
background light source and hence probe only the early phases 
of the starburst. Observations using a background source unre¬ 
lated to the host, such as a quasar or GRB, probe winds of all 
ages, which may have much different structure. The analysis 
of Theuns, Mo, & Schaye (2001) implies that winds at z ~ 3 
must fragment with Ca < 0.01 in order not to overproduce Lya 
forest lines. It is, however, unclear whether fragmentation of 
this degree will occur at the high redshifts we study given the 
smaller cosmic times available to the winds. 

Unfortunately, uncovering the properties of old winds will 
be difficult. When the host galaxy is used as the background 
source, the absorption lines themselves can constrain the struc¬ 
ture of the wind because saturated lines fix the fraction of the 
galaxy hidden by the absorbing medium (e.g., Heckman et al. 
2000). However, our background sources are pointlike and can¬ 
not be used in this way. One strategy to study the outflow struc¬ 
ture is to compare the strengths of the low and high ionization 
components of individual winds. A better strategy is to com¬ 
pare the absorption along neighboring lines of sight to close 
sources on the sky, such as quasar pairs (e.g., Bechtold et al. 





14 


1994; Crotts et al. 1994) or multiple images of a single gravita¬ 
tionally lensed source (e.g., Rauch et al. 2002). Lensed quasar 
images are particularly interesting because the lensing proba¬ 
bility increases rapidly with quasar redshift (Wyithe & Loeb 
2002a,b). Such studies can determine the projected extent of 
individual wind systems and the distribution of metal clouds 
within them. 

Our model has two important shortcomings. First, we ignore 
all effects arising from large-scale structure. If star-forming 
galaxies are highly clustered (as appears likely at very high 
redshifts), then their winds will overlap and decrease the fill¬ 
ing fraction of metal-enriched material. Clustered sources will 
also reside in overdense regions of the IGM, further reducing 
the wind bubble sizes. Qualitatively, both these effects will re¬ 
duce the number of intersected lines but increase their average 
strength. Clustering will also increase the variance in the num¬ 
ber of absorption lines per unit redshift, because those lines 
of sight that pass through overdense regions will have greatly 
increased numbers of intersections. Finally, the filamentary 
nature of large-scale structure will affect the geometry of the 
winds. Those winds expanding into filaments will stall earlier, 
while those expanding into voids will grow more quickly, re¬ 
sulting in anisotropic growth (Theuns et al. 2002). Thus most 
of the enriched volume will be inside voids, which would tend 
to increase the filling fraction. For example, with the extreme 
assumption that all winds expand into a medium with pg equal 
to one-tenth the mean density of the universe, the filling factor 
of enriched material increases by a factor ~ 3-5. On the other 
hand, if wind sources are embedded in overdense regions, the 
filling factor decreases even more dramatically because M^ax 
decreases rapidly. Unfortunately, a detailed treatment of clus¬ 
tering and the anisotropic expansion around filaments requires 
considerably more machinery than is developed here and is best 
addressed through numerical simulations. We have chosen a 
similar wind model to Springel & Hernquist (2002), so our re¬ 
sults should help to interpret the effects of the wind model pa¬ 
rameterization in such simulations. 

The second shortcoming of our model is the long-term fate of 
the winds. They are likely to fragment due to cooling instabil¬ 
ities or inhomogeneities in the ambient medium. If turbulence 
is present, the fragments will then mix with the surrounding 
IGM. Provided that this mixing occurs at approximately the 
sound speed, our estimates of the total filled volume will not 
be grossly changed; however, fragmentation could strongly af¬ 
fect the covering fraction of the metals (see above) and their 
velocity widths. 

There are at least two types of astronomical objects that can 
function as background sources. The most obvious are quasars. 
Bright quasars have been observed to z = 6.3 (Fan et al. 2002) 
and likely extend to z ~ 10 (Haiman & Loeb 2001). How¬ 
ever, the characteristic quasar luminosity is expected to decline 
sharply with increasing redshift (Wyithe & Loeb 2002c) and it 
would be difficult to find luminous quasars at much higher red- 
shifts. Quasars also suffer from complicated intrinsic absorp¬ 
tion, and isolating weak IGM metal lines may be challenging. 
Another possibility is to use gamma-ray bursts (GRBs) as the 
background sources. If GRBs are associated with massive star 
formation, some likely occur at extremely high redshifts (Lamb 

® See http://www2.keck.hawaii.edu:3636/realpublic/inst/nirspec/nirspec.html. 

See http://www.eso.org/instruments/isaac/. 

See http://www.stsci.edu/ngst/. 

See http://www.ngst.stsci.edu/nms/index.html. 


& Reichart 2000; Bromm & Loeb 2002). Cosmic time dilation 
and a favorable /T-correction imply that, at a given time after the 
burst in the observer frame, the decrease in brightness with z is 
mild. Ciardi & Loeb (2000) show that the 2/im flux of a GRB at 
z > 10 declines from ~ 200 pJy 1 hr after the burst to ~ 1 /rJy 
10 days after the burst. GRBs are also attractive sources be¬ 
cause of their intrinsic power-law spectra, from which it will be 
straightforward to find absorption lines. Furthermore, if GRBs 
occur in actively star-forming galaxies that drive winds into the 
surrounding IGM, the absorption spectrum can probe the host 
galaxy wind in detail. Because GRBs fade, they also serve as 
useful markers of faint galaxies at high redshifts. Subsequent 
deep imaging of the host could reveal details that complement 
the information from absorption line studies. 

Observing IGM metal lines at high redshifts presents a vari¬ 
ety of challenges. Most importantly, it requires relatively high 
resolution spectroscopy of faint sources in the near to mid- 
infrared (1/rm < A < 5pm). Resolving individual winds re¬ 
quires a spectral resolution R^c /Vexp ~ 6000 (50 km s“' /Vexp)- 
Such resolution is currently achievable with the NIRSPEC cam¬ 
era on the Keck telescope® and with the ISAAC camera on the 
Very Large Telescope.'® However, any ground-based obser¬ 
vation in this wavelength range will suffer from sky contami¬ 
nation that worsens significantly as the wavelength increases; 
thus, ground-based measurements will be limited to z < 10, 
and even then will require bright sources. The James Webb 
Space Telescope," due for launch in the next decade, will have 
the sensitivity to probe much fainter background sources with 
a spectral resolution R = 5000. The Mission Simulator^^ esti¬ 
mates that a source with /v ~ 1 pjy will require a ~ 16-24 hr 
integration to achieve S/N =10 per resolution element at wave¬ 
lengths \.5^pm. Thus spectroscopy of high-z GRB afterglows 
with this instrument will allow us to probe IGM metal lines at 
z > 10 in detail. 

Another challenge lies in identifying the lines. In the case of 
C IV and Si IV, the ratio of line strengths in the doublets will 
help to confirm the identity of the transition and hence the red¬ 
shift of the line. However, the low ionization states lack doublet 
transitions in the relevant wavelength range. One would there¬ 
fore need to identify multiple lines at each redshift in order to 
confirm the presence of an individual absorbing wind. The fea¬ 
sibility of this technique is then limited not by the strongest 
lines but by the weaker transitions. 

Finally, we note that absorption line observations such as 
those we propose can be combined with imaging surveys to 
learn even more about the enrichment process. Adelberger et al. 
(2002) have already shown the usefulness of such an approach 
for understanding the influence Lyman-break galaxies have on 
their surroundings. Though difficult, similar observations over 
a large redshift range have the promise to map out the history 
of galaxy feedback on the IGM. 
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Table 1 

Supernova Yields and Important Transitions 


Element 

Yx (Scalo) 

Yx (VMS) 

Ionization State 

(^) 

fosc 

C 

0.1 Mg 

4.1 Mq 

CII 

1334.5 

0.1278 




CIV 

1548.2“ 

0.1908 




CIV 

1550.8“ 

0.09522 

O 

0.5 

44 

Ol 

1302.2 

0.04887 

Si 

0.06 

16 

Si II 

1304.4 

0.094 




Si IV 

1393.8“ 

0.514 




Si IV 

1402.8“ 

0.2553 

Fe 

0.07 

6.4 

Fell 

1608.5 

0.058 




Fell 

2344.2 

0.114 




Fell 

2382.8 

0.300 


^Member of doublet 


Table 2 

Summary of Model Parameters 


Symbol 

Standard Value 

Definition 

f* 

0.1 

Star formation efficiency 

/esc 

0.25 

Fraction of supernova energy available for the wind 

/sw 

2/* 

Fraction of host gas entrained in outflow 

fd 

0 

Fraction of swept-up Hubble flow energy available to power wind 

fm 

0.1 

Fraction of swept-up material leaking into bubble interior 

I^SN 

126 Mq 

Mass of stars formed in order to produce one supernova 

£51 

1 

Explosive energy per supernova, in units of 10^' erg 

Yx 

See Table 1 

Yield of element x per supernova 

Ca 

1 

Areal covering fraction of metals 

Cm 

1 

Fraction of metals contained in given ionization state 




